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It has been believed since the paper by Fischler, Kogut and Susskindjl| that in 
QED2 at finite charge density the chiral condensate exhibits a spatially inhomoge- 
neous, oscillating behaviour. In this paper we demonstrate that this inhomogeneity 
is due to unphysical explicit breaking of the translational invariance by a uniform 
, background charge density. Moreover, we investigate in the context of a simple sta- 

^ \ tistical model what happens if the neutralizing background is composed instead of 

Q^' heavy, but dynamical, particles. We find that in contrast to the standard picture 



(D ' of the chiral condensate will not exhibit coherent oscillations on large distance 

scales, unless the heavy neutralizing particles themselves form a crystal and the 
density is high. 



X 



^ , PACS numbers: 

O 

a^ 

,Q ; I. INTRODUCTION 

o 

Over the past years there has been a lot of interest in the phase diagram of QCD at finite 
temperature and baryon density. The phase diagram would provide one with the answer 
Q-i! to the childish question, "What happens to matter when you heat it up or squeeze it?" 
^ I This question is relevant for the analysis of such extreme natural environments as the early 

universe or the dense interiors of neutron stars. The interest in the phase diagram of QCD 
has in turn sparked studies of models of QCD in non-trivial environments. 

Two dimensional Quantum Electrodynamics, QED2, commonly referred to as the 
I Schwinger model, has served as a playground for QCD theorists for many years. In the 
massless limit this model is exactly solvable and displays many features similar to QCD, 
most notably the generation of mass gap and the appearance of chiral condensate. The 
chiral condensate in this case is a manifestation of explicit chiral symmetry breaking by the 
axial anomaly and is generated in sectors with topological charge ±1. 

The first study of the Schwinger model at finite density has been performed in fT]. It 
must be noted that in this case we are talking about electric charge density, which unlike 
the baryon number density, is the zeroth component of a current associated with a gauged 
rather than global symmetry. Thus, to study 1-flavour Schwinger model at finite density 
one needs to introduce background charge in order to satisfy the Gauss law. The natural 
choice for such a background that was adopted in ^ is just a finite external charge uniformly 
smeared along the spatial direction. 

The massless Schwinger model remains exactly solvable at finite density. One of its most 
surprising features is that once an arbitrary small charge density is introduced, the chiral 
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condensate is no longer spatially uniform, but instead supports a plane wave structure [l| 



where p is the number density, is the topological angle of QED2 and {ipip)o denotes the 
chiral condensate at zero temperature, density and 6 parameter. Thus, the chiral condensate 
experiences oscillations with the period given by inverse number density. On the other hand, 
as long as the quark mass is vanishing the fermion density itself is uniform. So the model 
does not develop a conventional crystal but rather a "chiral crystal" . 

Once a finite quark mass m is introduced, the non-uniformity of the chiral condensate is 
translated into a non- uniformity of fermion density and to leading order in m/ejl|]. 



where ui = e/y^ and e is the gauge coupling constant. 

In this paper we will argue that the above conventional picture of oscillating chiral con- 
densate is unphysical, being a consequence of the introduction of an external static, 
neutralizing charge density. 

There have been a number of studies of the Schwinger model at finite density both in 
the Hamiltonianj^, 0] and path-integral formalism|3] since the work all of which have 
confirmed the oscillations of the chiral condensate. However, the ultimate reason for the 
breaking of the translational invariance in this model remains somewhat unclear. Indeed, one 
generally cannot spontaneously break continuous symmetries in 1 + 1 dimensions. Moreover, 
the ground state of the Schwinger model at finite density is unique (once the topological 9 
angle is fixed) and the lowest lying excitations are separated by a mass gap u. 

Yet, it is important to understand the precise reason for this phenomenon not only 
because it is curious by itself, but also as similar behaviour of the chiral condensate occurs 
in a large number of other models. In particular, both the 2d Gross-Neveau model and 
the two dimensional QCD in the large limit are believed to exhibit exactly the same 
spatial oscillations of the chiral condensate at finite density]^. The period of oscillations 
is again given precisely by the inverse fermion density. The stability of the chiral crystal 
against quantum fluctuations in these models is argued on the basis of the large N^. limit: 
once Nc = 00 one is allowed to circumvent the Mermin- Wagner theorem. Moreover, four- 
dimensional dense QCD in the large Nc limit is also expected to support a periodically 
modulated chiral condensate fi']. Thus, it would be valuable to first fully understand the 
origin of translational symmetry breaking in QED2, which is considerably simpler than the 
above zoo of models. 

We demonstrate, that the reason for this phenomenon is the presence of the background 
charge density, which leads to the inability to simultaneously maintain invariance under 
translational and large gauge transformations. Alternatively, in the path integral language, 
translational invariance is violated by sectors with a finite topological charge. These findings 
naturally explain the particular form of the chiral condensate at finite density and provide 
a more conclusive explanation for the loss of translational invariance than those present 
in the literature. We show that although both the chiral and translational symmetries are 
explicitly broken at finite density, in the massless limit a linear combination of them remains 
intact, which implies. 
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where 0{x) is an arbitrary local operator and q is the axial charge of O. 

Having thoroughly understood the reason why the uniform background density leads 
to explicit breaking of translational symmetry we ask the following question. Should we 
consider such a theory completely sick? More precisely, does the theory with the uniform 
background charge ever correctly model a theory where the neutralizing charge is heavy, 
but dynamical. Clearly, any theory where the neutralization of charge is performed solely 
by dynamical fields will not exhibit explicit breaking of translational invariance. Moreover, 
in the absence of "special arrangements" such as Nc = oo the translational invariance will 
not be broken in two dimensions spontaneously either. However, relics of the chiral crystal 
might remain intact on some finite, but large, distance scale. 

To answer the above question we consider a system where the neutralizing charge is 
modeled by dynamical classical particles of integer charge. We expect that this model 
corresponds to a theory where one fermion species is massless and the other is very heavy 
(of mass M), in the regime T <C M, e <^ M, where T is the temperature, provided that 
T is sufficiently large that the quantum effects for the heavy particles can be neglected. 
We integrate out the light degrees of freedom (photons and massless fermions) and are left 
with a classical statistical mechanics model for the heavy degrees of freedom. These heavy 
degrees of freedom have a size of roughly l/ut, interact via a Yukawa potential and should 
probably be identified with B-like mesons, consisting of one light and one heavy quark. 

We find that the chiral condensate in this model will not reproduce the standard picture 
of Q (see eq. (^), which exhibits for arbitrary density spatial oscillations with a density 
independent amplitude. Instead, the form of the condensate will depend crucially on the 
density and on the emergent dynamics of the mesons. In the regime where the model 
is tractable (i.e. when the mesons form a weakly interacting gas), the chiral condensate 
does not reproduce any of the features of eq. (H)). Instead, in the dilute limit p -C o;, 
the chiral condensate is uniform and its magnitude decreases slightly with density. The 
correlator {'ip^-Y-^{x)ip^-Y-4'iO)) does not experience any oscillations. In the high density 
limit, p ^ uj the chiral condensate decreases exponentially with density. The correlator 
{ip^^^{x)'ijj^-^'ijj{0)) exhibits oscillations with period on distance scales x <^ lu~^, 
which, however, disappear for x 3> uj~^. These oscillations on short distances are the only 
visible remnants of the chiral crystal in the gaseous regime. 

Thus, we shall argue that the chiral condensate has a chance to reproduce the plane wave 
behaviour (Q) on sufficiently large distance scales only if the density p ^ u and the heavy 
degrees of freedom themselves are close to crystallization. Unless these specific conditions 
are met, the uniform background charge approximation is inapplicable and the form of the 
chiral condensate ((H) found in 0, Q, 13 unphysical. 

II. WHAT'S NON-UNIFORM IN A UNIFORM BACKGROUND DENSITY 

This section is devoted to a detailed analysis of the reason for the appearance of the chiral 
crystal in a model with a uniform background density. The literature on this subject gener- 
ally supports the following argument present in the original paper If the spatial manifold 
is an infinite line R, one prefers not to introduce a background charge distribution that 
stretches across whole of M to avoid infra-red difficulties. Instead, one chooses a background 
charge density to be uniform in a certain finite region of the real line (say —L < x < L) 
and zero everywhere else. Once all the calculations are done one takes L ^ oo. Then the 
"small" explicit breaking of translational symmetry present in the form of the endpoints of 
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the charge distribution is carried by the long-range Coulomb forces across the whole system 
and leads to the chiral crystal structure (jT]). 

In principal, the above invocation of the long-range forces allows one to circumvent the 
general theorems on the lack of spontaneous symmetry breaking in 1 -|- 1 dimensions. How- 
ever, the above argument can no longer be directly applied once the spatial manifold is 
compactified to a circle with the background charge uniformly smeared along its length, 
apparently removing the "endpoints" of the charge distribution. We adopt precisely such a 
compactification of the spatial coordinate in what follows. 

Moreover, let us compare the situation to Schwinger model at zero density, where one 
observes breaking of the chiral symmetry. The modern philosophy regarding the origin of 
this phenomenon is that chiral symmetry is locally explicitly broken by the axial anomaly. 
Globally, one cannot simultaneously maintain invariance of the theory under chiral and large 
gauge transformations. Translating the last statement into the path integral formalism: axial 
charge is not conserved in non-trivial topological sectors. 

We now show that a similar picture holds for the breaking of translational invariance in 
Schwinger model at finite density. 



A. Hamiltonian Formalism 

We start with the Lagrangian for QEDi+i, 

L = -\f^.F^^^ + i^i^D^i, - m^^, D^ = d,- leA, (4) 

Local f/(l) symmetry of the theory takes the form, 

U{1) : i,{x) ^ e*"("V(a^), A^{x) A^{x) + ^d^a{x) (5) 

For the moment we work in Minkowski space, with the conventions 7^ = 7^7^, = 1. 
For definiteness, we take the spatial manifold to be a circle of length Li and pick the gauge 
where all fields obey periodic boundary conditions on this circle. 
The energy momentum tensor for this theory is, 

T'^^ = ^i^f^D^'i) - F>'^F\ - g^'^L (6) 

We have not symmetrized T^'^ as it is not essential for our purposes. 

Now let us couple the theory to a conserved external current iexti.^)i such that 9^jgxt = 0- 
The Lagrangian becomes, 

= L + i'^xt^fj- (7) 

Once this term is added, the energy momentum tensor satisfies, 

Q^j,,u ^pux^^^^^ (8) 

Clearly, an external current violates the conservation of energy and momentum. Now, let 
us take to represent a uniform, neutralizing charge density. 



iext = (-ep,0) 



(9) 
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so that p = N/Li, where is the total dynamical charge. 
The equation (jHJ becomes, 

df.Tf'^ = (10) 
d^T'^' = -epF (11) 

where F = Fqi, is the electric field. Thus, the uniform background charge density explicitly 
breaks spatial but not temporal translational invariance. In particular, defining the total 
momentum operator, 

P = jdx^T''^ (12) 

we obtain, 

—P = -ep / dx^F (13) 
If we integrate equation (|T!^ over time, 

AP^-ep/.^.F (14) 

We recognize the integral on the righthand side of eq. (fT^ as the topological charge of the 
theory. 

Thus, translational invariance is broken both locally and globally. One could try to 
redefine the energy momentum tensor as, 

f^"^ = T^'^ + epe^^A^ (15) 

and likewise the total momentum 

P = P + ep j dx^Ai (16) 

so that 

d,f^' = 0, = (17) 

However, the local current T^^ is clearly not a gauge invariant operator. The global object 
P is invariant under "small" gauge transformations characterized by a{Li) = a{0), where 
a{x) is the transformation parameter of eq. (0). However, P is not invariant under large 
gauge transformations U, 

U^{x^)U^ = e2"'"'/^i7/'(a;i), UAiU^ = A^ + ^ (18) 

whereby 

UPU^ = P + 2np (19) 

Thus, at finite background charge density, we cannot simultaneously preserve the invari- 
ance of the theory under both translational and large gauge transformations. The usual 
procedure, at least at zero density, is to formulate the theory in a way, which preserves the 
latter symmetry and to constraint oneself to states in the Hilbert space satisfying, 

U\9) = (20) 
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Then a finite translation with the operator P will take us out of the gauge invariant Hilbert 
space and into a different ^-vacuum: 

e'^y) = \e + 2npa) (21) 

Thus, for any local operator 0{x), 

{0{x))e = (O(0))e+2.px (22) 

This interplay between the 6 angle and the loss of translational invariance is clear from the 
expressions for the chiral condensate and baryon density (III),®- We would like to point 
out that we have not anywhere used the fact that our dynamical matter is fermionic. Thus, 
eq. ()22|1 would remain valid in a theory with any dynamical matter fields neutralized by a 
uniform background charge density. 

Note that a lattice subgroup of the translational group remains unbroken. Indeed, the 
operator, 

S = e'^f' (23) 

is invariant under the transformation (fTH)) . But S is an operator that performs translations 
by a distance a = - the average charge spacing. Thus, our theory will respect this 
symmetry as can be explicitly seen from (Q),®- 

The above discussion is precisely analogous to the philosophy behind the breaking of axial 
symmetry in Q-E-Di+i. Recall that the gauge invariant axial current, j^^ = ip'y'^'y^ip suffers 
from an anomaly, 

= -T^e^'i^/.. + 2tmilj^'^ix) (24) 

ZTT 

Equation is an operator identity and should not be affected by infra-red effects such as 
temperature or finite density. Let us define the following current, 

= T^^ - Txpf^ (25) 

Observe, that is a gauge invariant operator, satisfying, 

df.l^' = -27iimpip-f^tp (26) 

So in the massless limit m = 0, 

9^/^ = (27) 

Thus, at finite density, both the axial and the translational symmetries are broken. However, 
in the massless limit, a linear combination of them remains intact. Defining the global charge, 

/I 
dxU^x) = P -npQ^, ^<5 = (28) 

The conservation of Q dictates the structure of all "non-uniformities" provided that the 
symmetry associated with the conservation of Q is not spontaneously broken. Consider an 
arbitrary local operator of axial charge q, 

[Q^O(x)] = gO(x) (29) 

Then, 

(0(x)) = (fi|e"^«"0(a;)e^'3«|fi) = {Qle'^'f^^' 0{x + a)e-'''f"'^' \Q) = e"^'^'"^(0(a; + a)) (30) 
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{0{x + a)) = e~^"ip''{0{x)) (31) 



In particular, for the fermion bilinear ^i^^?/), q = =f2 and, 

{i^^-^i^{x)) = e^'-^^^^P^^m) (32) 

Thus, we see that the plane wave behaviour of the chiral condensate follows immediately 
from the structure of the theory. On the other hand, the density operator ipj^ip has q = 
and, therefore, does not display any non-uniformity in the massless limit. Thus, the equation 
(jnH) is in agreement with the explicit calculations of 

Once the quark mass m is non-vanishing the conservation of the current is explicitly 
broken (j26|) . Therefore, averages of local operators no longer need to satisfy the formula 
()31|). For instance, the fermion density (ip'j'^'ipix)) becomes non-uniform as can be seen from 
eq. 0. 

Before we conclude this section, we note that beside Schwinger model, both two dimen- 
sional chiral Gross-Neveu model and QCD2 are believed to display the structure in 
the large N limit In these theories both axial and translational symmetries are spon- 
taneously broken, but the linear combination (j28j) remains preserved by the ground state. 
Thus, the resulting picture is the same as in Schwinger model, but the formal reason for 
the appearance of a chiral crystal is very different. In Schwinger model, as we have shown, 
translational and axial symmetries are broken explicitly (by background charge density and 
by chiral anomaly). In the Gross-Neveau model and QCD2 these symmetries are broken 
spontaneously, with the theorems on the absence of spontaneous symmetry breaking in two 
dimensions circumvented due to = cxd. 



B. Path-Integral Formalism 

It is instructive to understand in parallel how translational symmetry breaking is realized 
in the path integral formalism. 
We go to Euclidean space with, 

Le = ^Ff^^Ff,^ + i'jf.D^i) + m^i) (33) 

In our notations 7172 = ^75 , ei2 = 1. 

We take the space-time to be a torus with < xi < Li, < X2 < 1^2. Physically, 
L2 = P = is the inverse temperature. Gauge fields on a torus fall apart into sectors 
classified by the topological charge. 



n = —Jd'xF (34) 

where F = F12. In a general topological sector, the gauge and fermion fields are not strictly 
periodic, but satisfy. 



V^(a;i, L2) = yi(xi)^(a;i, 0), A^{xi, L2) = A^{xu 0) - -9^Vi(xi)\//(xi) (35) 
V^(Li, X2) = ^2(^2)^^(0, X2), A^{Li, X2) = A^{0, X2) - -d^V2ix2)V^{x2) (36) 
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with Vi, V2 satisfying the consistency conditions, 

^(0)1^2(^:2) = V2{0)V,{L,) (37) 

The transition functions Vi, V2 in turn determine the topological charge n. 

For each n, we have some gauge freedom in choosing Vi,V2. For instance, one choice is 
to have fermions anti-periodic in the temporal (^2) direction, so that, 

= -1, V2{x2) = e2™^2/^2 (38) 

Let us recall that an external heavy static particle is inserted into the theory in the form 
of a temporal Wilson loop. For instance, the partition function in the background of m 
static charges located at points {xi} and with charges {epi}, pi G Z is, 

/m 
VAV^/jV^ Yl W{xi,pi) e-'^e™^ (39) 
i=i 

with 

W{x,p) = (^-Vi{xi)~'e''^^''^''^^''^y (40) 

We have inserted the prefactor (—1)^, so that in a gauge where Vi{xi) = —1, W{x,p) reduces 
to the standard form, 

W{x,p) = exp{ipe J A2{x,T)dT) (41) 

Expression (pn|) is completely gauge invariant and geometrically is the transport with respect 
to A along a temporal cycle, taken in representation p of the U{1) group. 

It is clear that once p ceases to be an integer the expression pUj) for W{x,p) becomes 
ambiguous (the only representations of the U{1) group are integral). This is not surprising - 
it is precisely for this reason that the existence of monopoles enforces quantization of electric 
charge^. In the present case the role of monopoles is played by 2d instantons. Similarly, 
it is problematic to generalize the prefactor in W involving the transition functions K to a 
continuous background charge distribution in a manifestly gauge invariant manner. 

We may still attempt to take the limit of a continuous charge distribution in a particular 
gauge. The choice Vi{xi) = —1 seems to be most suited for this purpose. As noted above, 
as long as we are working with integral charges in this gauge, the transition functions drop 
out of expression (00)). Now we can take the "continuum" limit, 

Z = j VAVipVip e^^/ '^'^Jr*^2 e-^e™'^ (42) 

where J2^*(a;i) is the static background charge density and (anti)periodic temporal boundary 
conditions on (fermions) gauge fields are assumed from here on. Expression ()42|) is not 
invariant under gauge transformations, which change these boundary conditions. 



However, the question of confinement of fractional charges in massless and massive Schwinger model has 
been discussed for agesj3,0,0|- This question has to be understood in the sense e'^^^i "^''''■^f' = ^^^pIo d- 
where D is the region of our manifold, such that dD = 7. The Wilson loop with the fractional charge is 
well defined only once we also choose D and is not independent of this choice. 
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Now let us address the question of translational symmetry breaking. First, to understand 
the root of the problem consider a fractional charge p situated at Xi = and move it across 
the artificial cut at Xi = ~ Li to Xi = 0+. Observe, 

W{Li,p) = e«ep/^2(il,T)a!r ^ ^iep J A2(0,T)dT J d2V2{T)V^ (t) ^ e^^'^^'^W {0, p) (43) 

Thus, the cut on the torus is visible to a fractional charge and invisible to an integer charge. 
Of course, there is nothing new in this result. However, it is precisely this fact that leads to 
translational symmetry breaking. 

Indeed, take an arbitrary local operator 0{x) and pick a > such that < x, x + a < Li. 
Let us compute the expectation value of (9(x) in the background of the charge distribution 
Then, 

(0(x + a)) = ^ J VAV^V^ 0{x + a)e''^ '''''' e-^e'"' (44) 

We make the following change of variables, 

ifj'ixi, X2) = ipixi + a, X2), A^(xi, X2) = Af,{xi + a, X2), < xi < Li - a (45) 
■?/''(xi,X2) = V2{x2)^{xi + a - Li,X2), Li - a < xi < Li (46) 

A^(xi, X2) = A^(xi + a - Li, X2) - ^df,V2{x2)V^{x2), Li - a < xi < Li (47) 

It is easy to see that ijj', A' obey the same boundary conditions as original variables ijj, A. 
Thus, 

{Oix + a)) = ^ J I?A'PV^'I)^'0'(x)e-2™^o'''='-U2^*(-i)e^^/'^'^'^r'^^ e'^e*"^ (48) 

where 

A'extf^ ^ _ j + a) < Xi < Li - a , . 

•^2 ^^^> - \ j|-*(xi + a - Li) Li - a < xi < Li ^^^^ 

is just the properly shifted background charge density. Thus, we get an extra factor, 

related to the amount of charge passing through the cut during our translation. As long as 
this charge is an integer, the factor (jSUj) is unity and the cut is invisible. 

Now suppose our background charge is uniformly smeared across the spatial circle, 
j|'^*(xi) = —p, where p = is the dynamical charge. The background charge density 
itself is invariant under shifts along the circle, jg^^* = and the factor becomes, 

(0(x + a))e = (0(x)) (51) 

in agreement with (j221)- Moreover, if we disentangle the contributions to (0(x)) coming 
from distinct topological sectors, 

(0(x))„ = e2™^^(O(0))„ (52) 

Thus, different topological sectors contribute to (0(x)) with different "harmonics." One 
effect of fl52p is that only the topologically trivial n = sector contributes to the partition 
function and the topological susceptibility vanishes even if the quark mass m is non-zero. 
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Finally, if the quark mass is vanishing then due to fermion zero modes, operators with 
axial charge q get a contribution only from sectors with n = — |. Thus, for m = 0, 

(0(x)) =e-"*^''^(O(0)) (53) 

in agreement with (j31|) . 

Before we conclude this section, we would like to note that in the massless limit ex- 
plicit calculations can be done for the theory defined by ()42|) on a finite torus. These are 
summarized in the appendix. In particular, 

(V^i±^^(x)) = ^e^^e2-^''^(^V')L„L. (54) 

where {iP4')li,L2 is the chiral condensate at zero density and 9 parameter on the torus of 
dimensions Li,L2. The result fj54|) generalizes the path-integral computation of j^, which 
was performed on an infinite Euclidean space to the case of a finite torus, where all the 
infra-red singularities are under complete control. Technically, the oscillating factor e^'^*^^ 
comes from averaging the fermion zero mode over torons (constant parts of field A^{x) = t^) 
in the presence of background charge. 



III. DYNAMICAL BACKGROUND 
A. General Remarks 

In the previous section we saw that a uniform background charge density explicitly breaks 
translational invariance. This effect may seem to be rather unphysical, in particular as it has 
to do with the Dirac string being visible to our background charge. But after all, the uniform 
background density is typically taken to model some heavy, but dynamical, particles. Once 
all fields are dynamical, it is clear that translational symmetry is not explicitly broken. 
However, one would like to ask whether any features of the chiral crystal discussed in the 
previous section remain. 

To answer the above question, we would like to analyze QED2 with two flavours. We 
take one fermion flavour to have charge e and vanishing mass and the other flavour to have 
charge ge, g G N and mass M ^ e. We want to analyze the problem at a flnite "isospin" 
density, with the heavy fermions neutralizing the light ones. We work at flnite temperature 
T. We will treat the problem in a "Born-Oppenheimer" like approximation. Namely, we 
flrst freeze the positions of heavy particles, treating them as static external charges, and 
integrate over the light fermions and gauge flelds. For instance, the partition function of the 
system in the background of N external charges situated at points {xi} is, 

Z{x,,..XN) = (Y[W{x„-q))i (55) 

i 

where the subscript / denotes integration over the light degrees of freedom. We then promote 
the external charges to dynamical degrees of freedom, treating them as classical particles. 
For example, the full partition function takes the form, 

n=0 ■ i 
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Here z = jc^*^^ is the activity, /i is the chemical potential and A = (j^)^ is the thermal 
wavelength. Similarly, the expectation value of some operator O involving light quark fields 
is, 

(O) = - ^ ^ / dxi..dxn{0 H W{x„ -q))i (57) 

n=0 ■ i 

We will often use the notation, 

(0){,^,,j = Z-i(xi,..xO(0 J]iy(x„g,))z (58) 

i 

We shall shortly see that after integration over the light degrees of freedom, the heavy 
fermions get dressed into meson like particles, consisting (in terms of quantum numbers) of 
q light and one heavy quark. So the effective theory ()56|) should be understood as describing 
classical dynamics of such mesons. We expect such an approximation to be valid as long as 
M 3> e, M T so that the heavy quark-antiquark pairs do not get excited either virtually 
or thermally. Moreover, we need T to be high enough that the meson gas/liquid is in a 
classical rather than quantum regime. In the dilute gas limit, we expect that the system 

2 

can be treated classically ioi T ^ 

As a first step to analyze the resulting system, we need to perform the integration over the 
light fermions and gauge fields. We shall work on a Euclidean torus of spatial and temporal 
lengths Li, L2 = T^^. The expectation value of a product of straight temporal Wilson loops 
in the massless Schwinger model has been computed in a number of works [8l.[lo|. The result 
is (see the appendix for a sketch of the calculation). 



i \ i.j 




(59) 



V{x) = ^y^T^e''^^ L.^oo ^gQ. 
Pi 

where u = e/y^ andpi = 27rm/Li, m & Z. Thus, our heavy particles (of like charge) inter- 
act via a two-body repulsive Yukawa potential with all three and higher particle interactions 
vanishing. It is also instructive to compute the charge density of light quarks, 

(^72^(a;)){x„q,} = -^Yl ^'^^^ ~ (^^) 

i 

It is clear from (jfilll that each heavy quark of charge —q is surrounded by a cloud of light 
quarks with a radius of roughly u!~^. The cloud has total charge q that screens the Coulomb 
potential of the heavy quark producing a meson, similar to the heavy-light mesons of QCD 
(such as the 5-meson). 

We will be most interested in the expectation value of the chiral condensate 
For static sources this is given by[l3] (we sketch the calculation in the appendix), 

(^^^^(^)){....} = U\i'i')L.M n (^(^ - ^^))"'' (62) 
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FIG. 1: a) The phase of the chiral condensate avgU{x) = aig{'il)^-^^ip{x)) in the background of 
a static unit charge placed at x = for = 0. b) The chiral condensate in units of (V'V')t in the 
background of two widely separated static unit charges for = 0. 



where 

U{x) = exp{2mV'{x)) ^'=°° exp(-7ris^n(x)e-'"l^l) (63) 

Thus, the introduction of static charges only affects the phase of the chiral condensate. 
Moreover, U{x) 1 for |x| ^ t^^^, so each static charge affects the chiral condensate only 
in a region of radius roughly uj~^ - the size of the meson. Notice that U{x) makes one loop 
on the unit circle in the complex plane as x winds around the spatial circle (see Fig. ^ 
a)). Thus, the phase of the condensate ()62|) winds by 27rN as x moves around the spatial 
circle, where N = — - qi is the total charge of the light fermions. So, the total winding 

number of {ip'^-Y~'^{^)) {xi,qi} is independent of the positions of the heavy quarks and, in fact, 
is the same as for the model with the uniform background charge density (j^lj) . However, 
the winding occurs in the vicinity of the heavy charges, over the radius of each meson, as 
opposed to the uniform background case, where the winding is uniformly smeared across the 
whole system.^ We expect this difference to be particularly important in the dilute limit 
p <^ uj when the distance between mesons is much larger than their size. In this regime, 
each meson keeps its individual features and {ip'ip){xi-q} ^ {ijjip)T cos{9) in the wide regions 
between the mesons (see Fig. Qb)). Here, {ipip)T is the infinite volume limit of the chiral 
condensate at zero density and 6 parameter and finite temperature T. Thus, the uniform 
background approximation is expected to fail badly in the dilute regime. 

It is instructive to see what happens to the chiral condensate if we arrange our heavy 
charges into a regular lattice, Xj = ja, qj = —q. Using (jU^ and taking Li oo, 

(V'^V^(x)){.,,,a = le^'i^^W e^'^(^) (64) 

where 0(x) is a periodic function with period x = a and, 

. sinh(u;(x - a/2)) 

= aTT^' < X < a (65 

smh(co'a/2) 



^ Technically, such a local nature of the result comes from a non-trivial cancelation between oscillating 
factors e^"'^^ originating from integration over different modes in the path integral (see appendix). 
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For c<ja ^ 1 we have a crystal of widely spaced mesons much like on Fig. ^h). In the 
high density limit, ua <^ 1, the screening clouds of heavy charges overlap and the individual 
mesons are washed out. Instead, we may approximate, 

~ 27rg - , 0<x <a (66) 

so that, 

ii^^-^^ix)) = \e'\-ire'-^^%i,i;)T (67) 

where p = q/a = N/L. Thus, in this limit we recover the uniform background approximation 
()54|) . However, note that the existence of coherent, long-range, oscillations of the chiral 
condensate ()67p is possible only if the dynamics governing the heavy charges are such that 
they organize a crystal-like state. For a one- dimensional statistical system interacting with 
a finite range potential (jUUI) a true crystal cannot form. However, the system may exhibit 
crystal order on some finite distance scale d ^ u!~^. In this case, we expect that the 
plane wave behaviour of the chiral condensate will also persist on the same distance scale d. 
Otherwise, if the mesons form a disordered, weakly-interacting gas the oscillations (|67|) will 
be washed out on distance scales x ^ as we shall show shortly. 

In fact, equation suggests that the external charges act as impurities, whose effect 
is to disorder the chiral condensate. If the impurities are in a weakly interacting regime, 
the disorder is "random". This is precisely the situation that we will analyze in the next 
section. 



B. Statistical Model 



Let us now make the heavy charges dynamical and analyze the statistical model (j^Bjl . 
Our main objective is to compute the chiral condensate 



n=0 i ij 

1 . - °° 1 /" 

= -e*^(^V^)T^— / dxi..dxnY\_iU{x - XiY - l)gn{xi, ..Xn) (68) 



where gn{xi, .-Xn) is the n-point correlation function, 

gn{xi, ..Xn) = z^— \ (ix„+i..rfx„+^exp(-- ^ g^eV(xi - Xj)) (69) 

m=0 ii=l 

Notice that the chiral condensate is sensitive only to short-distance properties of the corre- 
lation functions (7„(xi..x„) as the range of U{x)'^ — 1 is roughly u~^. 

We would like to perform the Meyer expansion in activity z. The leading term in the 
equation of state, as always, is, 

f3P = z* (70) 



14 



where P is the pressure and z* = z exp{—^Pq'^e'^V{0)) {z* includes the self-interaction energy 
of each meson). Then, 

= z^^{(3P) = z* (71) 
where is the density of heavy particles, 

QPh = P (72) 

As the range of the potential V is uj~^ and strength e'^V ~ u, the corrections to the equation 
of state for T ~ cij are suppressed in the Meyer cluster expansion by powers of z*/uj. So 
for, p <^ uj, T uj our system behaves like a weakly-interacting gas of mesons. Moreover, 
in this regime, at leading order the n point function on distances x ~ uj~^ scales as 2;*" so 
that terms in (jUHj) involving gn{xi, --Xn) are suppressed by {z* /uj)"'. The leading correction 
to the chiral condensate comes from the n = 1 term. Recalling, gi{x) = ph, 

(^^^^) ^ ^e^'(^^)T(l + Pk f dx{U{xy - 1)) (73) 



Performing the integral. 



where, 



(^^) ~ (1 _ P!Lu{q)){i,^)^cos{e) (74) 

UJ 



u{q) = 2 / dt^- -AiZ = 2(log(7rg) + 7 - «(7rg)) (75) 

Jo ^ 

Thus, we have calculated the first correction in density to the chiral condensate in the regime 
of a weakly interacting dilute meson gas. The result ()74|1 is reminiscent of the behaviour 
of the chiral condensate in A'^c = 2 QCD at small baryon densityjllj, in A^^ = 3 QCD at 
small isospin density l2| and in the "dilute" nuclear matter [isf. In all of these theories, one 
thinks of the system as being composed of a dilute gas of particles M (diquarks for Nc = 2 
QCD, pions for Ac = 3 QCD, nuclons for nuclear matter) and obtains, 

iij^), ^ (V^V^)o + p{M\^^\M) (76) 

where \M) is a one particle state in vacuum (with normalization {M{p)\M{p')) = (27r)'^(5'^(p— 
p')). So, we identify, 

{q\^i^\q)T,e = [ dx{U{xy - 1) {ijij)TCOs{e) = --M(g)(V^V^)rcos(^) (77) 
J ^ 

where |g) denotes our heavy-light meson state. 

So far we have concentrated on the region T uj where the criterion for the applicability 
of Meyer's expansion was p <^ u. Now, let's analyze the low temperature regime T ^ cu. 
In this case, the interaction effectively becomes hardcore of range uj~^log{uj/T), so for the 
Meyer expansion to be vahd, we need z* <C u/ \og{u/T). If this condition is satisfied, the 
chiral condensate at leading order in p is again given by ()74|) . Moreover, we actually expect 
the expressions (f7H) . (f77j) to remain valid in a wider range p ^ oj down to the extreme 
quantum regime at T = 0, based on general phase-space arguments. 
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Finally, let's study the high-temperature regime T ^ uj. In this case, it can be shown 
that the corrections to the ideal gas equation of state (f7n|) are suppressed by powers oi z/T. 
In particular, if a; -C p -C T we are still in the weakly interacting (but not dilute!) regime. 
In this case it is convenient to rewrite (jUS|) as, 

(V' ^ \^ = ^e'^{i>^lj)TeW \y ]-T I dXi-.dXn I \ {U{x - XiY - l)gn{Xi, ..Xn)c 



2 



\n=l ' i 



(78) 

where gn{xi, ..Xn)conn denotes the fully connected n-point correlation function. It is easy to 
show that the terms in the exponent involving the n-point correlation function are suppressed 
by {z/T)^~^ compared to the leading term and, 

(i/jiP) ^ exp (V^V^)t cos(^) (79) 

The expression fl79|) agrees with (ff^ in the dilute limit p <^ uj. In the dense gas limit, 
uj <^ p , the chiral condensate exponentially decreases with density. Note that for T ^ a; 
the chiral condensate at zero density is already exponentially suppressed with temperature 
compared to T = (see eq. ()A28|) ). 



C. Correlation Functions 

To answer the question of whether any remnants of the oscillating behaviour ()54|) exist 
in our model, the computation of the chiral condensate presented in the previous section is 
not sufficient. Indeed, translational invariance implies that the chiral condensate is uniform. 
Instead, we must compute the static correlation function, (5+ (x)S'_ (?/)), S±{x) = ilj^^^tlj{x). 
We would like to see on what distance scales this correlation function exhibits plane wave 
structure (j33)). Integrating out light degrees of freedom, 

(S+(x)5_(y)){.^,,,} = S{x - y) \{{U{x - x,))-'^^{U{y - x,)^ (80) 

i 

where, 

S(x-y) = (S+(i)S_(!/)>L„fa = i|W'>l.,LjV''°-'l-'> (81) 

G"W = 7VE^»*''- P=f^.^)."..".eZ (82) 
L1L2 ^ +oj^ \ Li L2 J 

Thus, 

{S+{x)S^{y)) = S{x-y)F{x-y) (83) 



F(x - ?/) = exp I ^ ^ I dxi..dxnW{U{ 

\ri=l ^' ^ i 



X ~ XiYU{y - Xi) -l)gn{Xl, ..Xn)connj (84) 

So the correlation function factorizes into two pieces. The first, S{x — y), is just the corre- 
lation function in vacuum. The second, F{x — y), contains the finite density information. 




As noted in the previous section, in the regime where the Meyer expansion is apphcable, 
we may truncate the series in the exponent of (j85|) at the leading [n = 1) term, 

F{x -y)= eP'^f^^-y^ (85) 

where, 

fix -y) = j dx,{U{x - x^YU{y - xi)-" - 1) (86) 
The function /(x) is plotted in Fig. |21 It is easy to see that for \x — 

f{x -y)^ j dxi{U{x - xiY - 1) + y dxi{U{y - xiY" - 1) = ~u{q) (87) 

Thus, {Sj^{x)S_{ii)) will not exhibit oscillations for \x — y\ ^ uj~^ . 
Note that S{x — y) ^ iK'^V')Tp for \x — y\ oo so, 

{S^{x)S^{y)) -\{i,i^)T? eM-'^-n{q)) = {S+){S-) (88) 

4 UJ 

and the correlation function clusters. 

We may also investigate the short distance behaviour. Expanding f[x) in a Taylor series 
in ux, 

f{x) = -{2'KiquJx - -■K'^{qujxf + ...) (89) 
UJ 2 

Hence, for \x\ <^ cj^^, 

2 

(5+(x)S_(0)) ^ exp(27ripx - '^^{ujxf)S{x) (90) 

The above equation clearly exhibits the plane wave behaviour with period x = p^^. However, 
recall that eq. (fnUjl . is vahd only for |x| <ti uj~^. Thus, in the dilute limit p <^ u, no full 
oscillations appear and, in fact, eqs. . are more appropriately written as, 

{S4x)S40)) = il + pHfix))Six) (91) 

In the dense gas regime, p ^ u, the oscillations are, indeed, present on short distance scales, 
however, as eq. shows, they become damped for x > {^)2uj~^ and disappear altogether 

for X ^ uj~^. Moreover, these oscillations modulate the zero-density correlator S{x), which 
itself has a quite non-trivial behaviour for distances x < u^^. 
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IV. CONCLUSION 

In this paper we have analyzed some puzzles related to the Schwinger model at finite 
density. We have shown that the well-known plane-wave behaviour of the chiral condensate 
is a consequence of explicit breaking of translational invariance by a background charge 
density. Similarly to the non-conservation of axial charge, the non-conservation of total 
momentum is globally saturated in sectors of non-trivial topological charge. In fact, the 
breaking of translational symmetry at finite density is a much simpler phenomenon than the 
breaking of chiral symmetry by the anomaly as the former appears already on the classical 
level, while the later is a purely quantum phenomenon. 

In the second part of this paper, we have explored the question: "What happens if the 
uniform background density is replaced by a dynamical, but heavy, field?" To answer this 
question, we have analyzed a statistical model in which the heavy neutralizing charge comes 
from an ensemble of classical particles. We have shown that the effect of heavy charges is 
to disorder the chiral condensate. In the regime where the gas of heavy charges is almost 
ideal, the chiral condensate is spatially uniform and decreasing with density. For the charge 
density p -C a; the "disorder" is weak and we compute the leading density correction to the 
chiral condensate ()74|) . In the dense gas regime p ^ uj the disorder leads to an exponential 
suppression of the chiral condensate. In both of these regimes, the condensate does not 
exhibit any oscillations on distance scales x 3> u!~^, as is clear from computing the correlator 
(S'+(x)S'_(0)) = {'^p^-Y-^{x)^jJ^-^^p{0)). The only remnant of oscillatory behaviour comes 
at high density p ^ u in the short distance behavior of {S+{x)S-{0)) for x <C 

In fact, we have argued that the only way for the oscillations to survive on distance scales 
larger than uj~^ is for the system to be in the high density regime p ^ uj and the heavy 
charges to crystalize. In the dilute regime, p -C a; we do not expect to recover the oscillatory 
behaviour even if the heavy charges were to crystalize. 

So clearly the uniform background approximation generically does not accurately model 
the situation where the neutralizing charge is dynamical, rendering the results of 0, 0, 
0| unphysical. Indeed, we expect such an approximation to work well if the light and 
heavy neutralizing charges are largely decoupled from each other (e.g. valence electrons in 
a metal). However, in the Schwinger model the light fermions are very strongly coupled to 
the neutralizing charges producing heavy-light mesons. The approximation fails particularly 
badly in the dilute phase, where the distance between the mesons is much larger than their 
size. To apply the uniform background charge approximation here would be akin to treating 
the nucleii in a dilute atomic g uniform. 

We conclude by noting that we certainly have not analyzed the entire phase-diagram of 
the two-flavour heavy-light Schwinger model. We have treated the gas of heavy-light mesons 
classically and have not touched upon the quantum regime at all. Neither have we analyzed 
the regime where the classical system is far from an ideal gas limit and the interactions 
between mesons are important. These regimes are subject to further investigation and, in 
fact, have a higher chance of exhibiting the plane- wave behaviour of chiral condensate on 
distance scales larger than uj~^ than the "random disorder" case considered here. 
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APPENDIX A: EXPLICIT CALCULATION OF THE CHIRAL CONDENSATE 

AT FINITE DENSITY 

The purpose of this appendix is to perform the calculation of the partition function and 
chiral condensate at finite background charge density on a Euclidean torus. In case when 
the background charge is in the form of discrete charges (Wilson loops) this computation 
has been performed beforej^^^. For a uniform background charge density, the calculation 
has been done on an infinite Euclidean plane Here, we keep the size of our torus finite 
throughout the calculation, gaining complete control of all the infra-red subtleties. For a 
detailed study of the Schwinger model on the torus (at zero density) see 3 • 



We work in a gauge where the (fermion) gauge fields are (anti) periodic in the temporal 
direction, with the transition functions 



Vl(xi) = -1, V2{X2) = e^™/^^ (Al) 

We decompose the gauge fields as, 

A^ = t^ + d^a + e^M + Al (A2) 

where a and h are both periodic fields on the torus orthogonal to unity (J d?x ex = j cPx h = 
0). The variable is the so-called toron field and plays a crucial part in all the calculations, 
is effectively an angular variable, with ~ + G Z. We shall consider only the 

case where the total background charge is integral, J dxij2^'^{xi) = —N, e Z so that the 
angular nature of is unspoiled. is the "instanton" field in the n-th topological sector. 
We choose, 

A-O. ^-H^ (A3) 

which obeys the periodicity conditions 

First, let's compute the partition function. 



with the normalization Z = 1 for jg^* = 0. For vanishing mass, only contributions from the 
trivial topological sector survive (recall that there are precisely |?t,| fermion zero modes in a 
sector with topological charge n with 75 = sgn{n)) and, 

Z= f Z^A det(7^D^)e-5/'^'-^'e'^-^'^'"^'2"^^ (A5) 



The (regularized) Dirac determinant is given by (see jl4| and references therein), 

det'(Li(7^D^)) = det Arexp(r(6) + r(t,n)) (A6) 
where det' denotes the determinant with the zero mode contributions deleted and, 

r(6) = f d?xd,hd,h (AT) 



2 







\n\ log 









(AS) 
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Here, 9 = 9i + 162 = -j- = iLz TV is the matrix of zero mode inner products, 



2lT ' 



n„ = j <e.-Ai^xM (A9) 

^^.{x) = e^^°(^')e^^^^(^)x°(a;) (AlO) 

where x^{^) ^ire the orthonormal zero modes of the operator Dq = '^f_i{d^ — ie{t^ + A^). 
Thus, 

/ dtidt2 e-»e^*2^2gr(i,0) J pfe ^ie^ d?xd^jf=\x^) b(x)^-S(h) 

where, 

s{b) = ^J d\b{x){-d^){-d^ + uj^)b{x) (A12) 

Performing the integral over the toron fields we obtain, 

-exp(-L2-— ) (A13) 



we recognize ei;' = as the Fermi energy of free massless Dirac fermions in 2d at fermion 
number N. 

The integration over b field gives, 

JVbFsW) 

with the propagator. 



exp Q y d^xd^yr]{x)G{x - y)r]{y)^ (A14) 



= i^(G'o(x)-G'^(x)) (A15) 
G.(x) = -i-E^^^^^ P=fe^ym„^,ez (A16) 
So we find. 



exp J dxdyjT\x){V{x-y)-^)jT\y)^ 

(A17) 

with given by (jHOI). Combining the global and local pieces ()A13|1 . ()A17|1 . 

Z = exp (^-L2^ j dxdyjT\x)V{x - y)jT\y)^ (A18) 

For a uniform background charge density, jg^j ~ contribution to ()A18|) comes 

only from the global piece and is given by ()A13|1 . On the other hand for discrete integral 
charges, jl^t{x) = Qi^ix - Xi) and, 

Z={Y[ W{x„ q,)) = exp ( "^2^ mje^V{xi - x,) J (A19) 

i \ ij / 
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Now, let's compute the chiral condensate, {ip^^^ip{x)). For m = 0, it receives a contri- 
bution only from the topological sector with n = 1, 

2 Jn=l 2 

(A20) 

where ^ is the normalized zero mode of the operator '~f^D^. We have, 

^{x) = Af-^xix) = Ar-5e*^"(^')e'=^5^(^)x°(a;) (A21) 

As noted, x° is the normalized zero mode of the Dirac operator in the background of toron 
and instanton fields, 

j,{d,-ie{t, + A;))x%x) = (A22) 
obeying the boundary conditions . . In the n = 1 sector we have a single zero mode, 

X'{x,, X,) = (^y^ ' e--l-l''?e^(''"^"^-"^)+'^"^V4(vr(xi + tx^) - 7r^r|r)e+ (A23) 

where Xj = Xi/Li and ^+ is the spinor with positive chirality 75^+ = 
Thus, 

J x>b e*^ / d-^^'dihix'^) b(x')^-s{b) 

where 5*0 = gilli^ is the "bare" instanton action. Performing the average over the toron 
fields, 

jdt,dt,L,xlx.{x)e-^'^^^e^^'''^ _ 1 ..n..^Nx.,l, 

/rftidt2e-^*^^^er(*.o) " T,"^ ^^^^"^^ " ^^^^^ 

Taking the average over the b field with the help of ()A14|) . 

fj)!^ 2ebix) iefd''x'dij2{x[)b{x') -Sib) . , . f P \ 

, L = 6^^^^°) exp 2ze^ / ^^^^'^(a; - x')d,jr\x[) = 

= e2^'^(°) exp (^2m J dx[ (^^^ - l-sgnix, - x[) - V'ix^ - x[)^ jf* K)) (A26) 
Now, combining eqs. ()A25|) . ()A26|) . 

i^^^Hx)) = le'\-lf{i^i^)L,M^^P f-2« / dx[ (l-sgn{x, - x[) + V'{x, - x[)] jr\x[) 



(A27) 

where the chiral condensate in the absence of external charge and at 6' = is given by. 



Li x^o ir\x\ I 



2 „-S'n„2/_N _2e2G('0^ ^ ^-2nG,.,(x) I ^^'^ ^ T — > 0, Li ^ OO 

2Te~~ T — i> oo, Li ^ oo 

(A28) 
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If the background charge density is uniform, 

In this case the oscillating factor comes solely from the integration over the toron fields 
()A25|1 . On the other hand, if the external charges are discrete and integral. 



i \ i 



(A30) 

and the long range oscillating factor is canceled between the global ()A25|) and local ()A26|) 
parts. 

A similar computation can be performed to obtain the result ()8()|1 for the correlation 
function of chiral densities. 
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